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ABSTRACT 

We present XMM RGS and Chandra LETG observations of the blazar, H 2356-309, located be- 
hind the Sculptor Wall, a large-scale galaxy structure expected to harbor high-density Warm-Hot 
Intergalactic Medium (WHIM) . Our simultaneous analysis of the RGS and LETG spectra yields a 3<r 
detection of the crucial redshifted O VII Ka line with a column density (> 10 16 cm -2 ) consistent with 
similar large-scale structures produced in cosmological simulations. This represents the first detec- 
tion of non-local WHIM from X-ray absorption studies where XMM and Chandra data are analyzed 
simultaneously and the absorber redshift is already known, thus providing robust evidence for the 
expected repository of the "missing baryons" . 

Subject headings: X-rays: diffuse background — X-rays: galaxies: clusters — cosmology: observations 
- cosmology: diffuse radiation — galaxies: BL Lacertae objects: individual : 
H 2356-309 



1. INTRODUCTION 

At most 50% of the total baryonic matter in the nearby 
universe can be accounted for by the amount of luminous 
baryons revealed by observations of stellar light, narrow 
Lya absorption systems, a nd X-ray emission fro m hot 
gas in galaxy clusters (e.g.. iFukugita et "aLl [l998). Cos- 
mological simulations predict that most of the "miss- 
ing baryons" (30-50%) reside in low-density plasma in 
large-scale filamentary structures betwee n galaxies (e.g., 
ICen fc Ostriker1ll999t [Dave et al.ll2001h . The location 
and predicted temperature range (10 5 — 10 7 K) of this 
plasma motivated its name as the "Warm-Hot Intergalac- 
tic Medium" (WHIM). 

Ultraviolet observations of O VI absorption lines in 
background quasar spect ra provided the first clear detec- 
tion of the WHIM (e.g., ISavage et all 119981 ; iTripp et all 
l2006l and references therein). The O VI lines probe 
lower temperature WHIM (ss 10 5 K), while most of the 
WHIM gas is expected to exist at higher te mperatures 
obser vable only at X-ray wavelengths (e.g., iDave et all 
2001). Important evidence for WHIM in X -ray emis- 
sion has been provided by, e.g., correlation studies of 
X-ra y emission and galaxy pos itions in superclusters 
(e.g. JZappacosta et alj2 002, 2005]) and by X -ray imaging 
studies of the highest density WHIM near massive clus- 
ters (e.g.. IWerner et af1l2008l and references therein). 

Because of the difficulty in treating the foreground X- 
ray emission from the Milky Way, it is important to 
corroborate these detections with X-ray absorption-line 
studies which are largely insensitive to this issue. Typi- 
cally, X-ray absorption studies focus on the O VII Ka line 
(A = 21.6 A), because it is expected to be the strongest 
line at X-ray energies owing to the anticipated WHIM 
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temperature and the relatively high cosmic abundance 
of oxygen. 

To date there is no detection of WHIM in X-ray ab- 
sorption that is gen erally accepted by the astronomi- 
cal community (e.g.. iRichter et aTll2008l and references 
there in). The case of M kn 421 illustrates the prob- 
lem. iNicastro et al.l (|2005[) first reported highly signif- 
icant detections of two absorption systems along the 
Mk n 421 sight lin e usin g a deep Chandra observation. 
But iKaastra et al.l (|2006f ) showed that the statisti c al sig - 
nificance of the lines reported by INicastro et al.l (2005) 
was greatly overestimated because the redshifts of the ab- 
sorbers were not known a priori. Furthermore, a deeper 
XMM observation of Mkn 421 did not confirm the pres- 
ence of the claimed absorption systems (|Rasmussen et all 
12001 . 

To address these issues, we consider a different obs erva- 
tional approa ch proposed by iKravtsov et ail (|2002f ) and 
pioneered by iFujimoto et ail ( 20041 ) in their absorption 
study of the Virgo cluster. Rather than performing a 
so-called "blind" search for previously unknown inter- 
vening WHIM structures in front of the brightest back- 
ground AGN, we perform a "targeted" search by choos- 
ing a region of sky that provides an optimal combination 
of known large-scale structure and a bright, intrinsically 
featureless, background source. Because the redshift of 
the foreground structure is known a priori, the statisti- 
cal significance of a detected line is much higher than for 
a blind search. By focusing on large-scale structure, we 
should observe relatively high density WHIM, therefore 
aiding a detection. We restrict our study to background 
blazars because of their intrinsically featureless spectra 
compared to other AGN. 

In this paper we present XMM and Chandra grating 
observations of the blazar H 2356-309 (z = 0.165) lo- 
cated behind the Sculptor Wall (Figure Q}, a southern 
super-str ucture filled with ma ny groups and clusters of 
galaxies (|da Costa et al.|[l~994f) . The part of the Sculptor 
Wall intercepted by the sight-line of H 2356-309 repre- 
sents a redshift range of 0.028-0.032. Although the obser- 
vations were designed as Targets-of-Opportunity to catch 
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Fig. 1. — Sky map (top) and wedge diagram (bottom), each ex- 
pressed in R.A., of the region of the Sculptor Wall where the blazar 
H 2356-309 is located. Galaxy redshifts are taken from NED. The 
(blue) dashed lines in the upper panel define the range of declina- 
tion of the galaxies displayed in the lower panel. Conversely, the 
(blue) dashed lines in the lower panel define the redshift range of 
the points in the upper panel. The blazar line-of-sight is indicated 
by the (red) dotted line in the wedge diagram. 

H 2356-309 in outburst, the apparent outbursts were not 
long enough to enable the follow-up (10-40 hours after 
triggering) to catch the blazar in a high state. The ob- 
servations were performed four months apart because of 
different triggering criteria used for the different satel- 
lites. 

Despite H 2356-309 being observed in its low state 
(~ 10~ n erg cm" 2 s _1 0.5-2.0 keV flux), we detected 
a candidate WHIM O VII resonance line in the Sculptor 
Wall from a simultaneous fit of the XMM and Chandra 
data. Although we examined other candidate lines from 
other parts of the spectra using all the XMM and Chan- 
dra data, here we present only the results for the O VII 
line because only for it do we achieve a detection of at 
least 3cr significance. 

The paper is organized as follows. In SJ21 we present 
the observations and describe the data preparation. The 
spectral fitting, modeling approach, and systematic er- 
rors are discussed in <JJ] We present our conclusions in 



2. OBSERVATIONS 

On June 2, 2007 XMM observed H 2356-309 (ObsID 
0504370701) for approximately 130 ks with the Reflec- 
tion Grating Spectrometer (RGS) as the primary instru- 



ment. The RGS is comprised of two nominally identical 
sets of gratings, the RGS1 and RGS2. But due to the 
failure of one of the CCD chips early in the mission, the 
RGS2 lacks sensitivity over the energy range (20— 24 A) 
relevant for the O VII lines. Consequently, we focus on 
the RGS1 in this paper. We analyzed the RGS1 data 
with the most up-to-date XMM-Newton Science Analy- 
sis Software (SAS 8.0.0) 6 along with the latest calibration 
files. 

We generated a light curve of the background events 
using CCD number 9; i.e., the CCD that is close to the 
optical axis and is most likely to be affected by the back- 
ground flares. Inspection of the light curve does reveal 
a flare near the end of the observation, which we re- 
moved, resulting in a clean exposure of 126 ks for the 
RGS1. Using these good time intervals, we reprocessed 
the RGS1 data to produce the data files required for spec- 
tral analysis; i.e., the response matrix, the background 
spectrum, and the file containing the spectrum of H 2356- 
309 (which also contains background). We restrict our 
analysis to the first-order spectra, because the second or- 
der does not cover the relevant O VII line energies and 
its count rate is much lower. 

Chandra observed H 2356-309 for 95.5 ks with the Low- 
Energy Transmission Grating (LETG) and the HRC- 
S detectors, which offer the best compromise between 
sensitivity and spectral resolution in the wavelength 
range of the O VII lines. We reduced the data using 
the standard Chandra Interactive Analysis of Observa- 
tions (CIAO) software (v4.0) and Chandra Calibration 
Database (CALDB, v3.4.0) and followed the standard 
Chandra data reduction threads . To ensure up-to-date 
calibration, we reprocessed the data from the "level 1" 
events file to create a new "level 2" file for analysis. We 
applied the lastest HRC gain map and pulse-height fil- 
ter for use with LETG data 8 . This procedure removed 
a sizable portion of background events with negligible 
X-ray event loss. From inspection of the light curves ex- 
tracted from source-free regions of the detector, we con- 
cluded that the observation was not affected significantly 
by background flares. 

We used the CIAO software to produce the files re- 
quired for spectral analysis; i.e., the response matrix, the 
background spectrum, and the file containing the spec- 
trum of H 2356-309 (which also contains background). 
Unlike the RGS, different orders cannot be separated 
from the LETG-HRC-S spectra, and so a given spec- 
trum contains all orders. Consequently, by default we 
use a response matrix that includes information up to the 
sixth order. In M3.5I we compare results using a response 
matrix that includes information only on the first-order 
spectrum. 

We re-binned the spectra to optimize detection of 
the absorption lines. After some experimentation, we 
achieved this by requiring a minimum of 75 and 40 counts 
per bin respectively in the RGS1 and LETG spectra. 

In Figure [5] we present the background-subtracted 
Chandra and XMM spectra. (The background subtrac- 
tion is achieved in the normal way in xspec via the " data" 
and "backgrnd" commands.) Two candidate absorption 

6 http: / /xmm.esac.esa.int/sas/8.0.0/ 

7 http:/ /cxc. harvard.edu/ciao/threads/ 

8 http:/ /cxc. harvard.edu/contrib/lctg/GainFilter/ 
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Fig. 2. — Background-subtracted XMM (black) and Chandra (red) spectra of the blazar H 2356-309 over the wavelength range 2f .0 — 22.5A 
used to study the candidate O VII lines. The spectral models represent a power-law continuum, two absorption lines broadened by the 
instrumental resolution and the Voigt function, and foreground Galactic absorption from cold gas. The wavelength positions of the two 
absorption lines, nearby and Sculptor, are indicated. 



lines in each detector are immediately apparent upon 
visual inspection. The line near 21. 6A corresponds to 
O VII for z = 0, consistent with nearby WHIM in the 
Local Group or hot gas in the Milky Way. However, the 
other candidate line near 22. 3 A corresponds to O VII for 
z w 0.03, consistent with the redshift of the unvirializcd 
structure in the Sculptor Wall. (The absorption line near 
2 1.8 A in the RGS1 is a known instrumental feature.) 

3. SPECTRAL FITTING 

We perform ed spectral fitting using XSPEC vl 1. 3.2a; 
( Arnaud 1996) and minimized the Cash C-statistic (|Cas^ 
I1979I ). rather than x 2 , to obtain unbiased parameter 
estimates from the P oisson-distributed spectral data 
([Humphrey et al.1 12008) . In all the fits the foreground 
absorbing Hydrogen column density from cold gas in the 
Milky Way was fixed at the value (1.33 x 10 20 cm -2 ; 
iDickev fc Lockmanl I1990D using the phabs model in 

XSPEC. 

3.1. Models 

A common practice in X-ray studies of the WHIM is 
to characterize the properties of a candidate absorption 
line by subtracting a gaussian component from a model of 
the continuum. The equivalent width is computed from 
this composite model, from which the column density of 
the absorber is inferred. This procedure is well suited 
to the case where the absorber is optically thin (i.e., lin- 
ear part of the curve of growth), which we show is not 
the case for the candidate Galactic and extra-galactic 
WHIM absorbers in H 2356-309. In these cases a com- 
posite continuum-minus-gaussian model takes negative 
values at the line centers. Consequently, we constructed 
a physical line absorption model as follows. 

Let I(E; 0) be the monochromatic specific intensity of 
photon energy E incident at coordinate s = assum- 



ing plane-parallel geometry. The emergent radiation at 
coordinate s is then, 



I(E;s)=I(E;0)e-^ E ' s \ 
where the optical depth 9 is, 

t(E;s)=N(s)—ME), 
m e c 



(1) 



(2) 



The quantity, N(s) = nds, is the column density of ab- 
sorbing atoms, and / is the absorption oscillator strength 
of the spectral line transition. The line-profile function, 



4>(E) = 



=H(a,u), 



AE D ^ 

is expressed in terms of the Doppler width, 

AEd = Eq~, 
c 

where Eq is the rest-frame energy of the line, and 



(3) 



(4) 



6=7^=32.24 



T 16 km 

s 



10 6 K A 



— , (5) 



is the Doppler 6-parameter. We have expressed b in terms 
of an oxygen atom of 16 atomic mass units. The quantity 
H(a, u) is the Voigt function which combines the effects of 
natural and Doppler broadening through the parameters, 



hT 



At: AE D ' 



AE 
AEd 



(6) 



9 We do not correct for stimulated emission because (1) the effect 
is small for the resonance O VII line for expected WHIM temper- 
atures, (2) the WHIM temperature is not precisely determined by 
our fits, and (3) the (1 — exp[— E/KgT]) correction term applies 
assuming a thermal plasma, which may not be strictly valid for the 
WHIM. 
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where, T is the Einstein-A coefficient of the line transi- 
tion, and AE = E — Eq. We use d the imple mentation 
of the Voigt function provided by IWellsl (|1999l ) 10 . We 
implemented the line absorption as a local multiplica- 
tive model in xspec, such that for a given energy bin 
of an input response matrix we evaluated the average 
of exp[— t(E; s)] over 10 equally spaced energies within 
the bin. (The variation of I(E; 0) over an energy bin in 
the RGS and LETG response matrices is negligible.) For 
the resonance O VII Ka line we used Eq = 0.574 keV 
(21. 6 A), / = 0.696, and T = 3.3 x 10 12 s^ 1 dVerner et alJ 
fl99"l . 

We find that a power-law model is a good description 
of the continuum over 21.0-22.5 A (0.5515-0.5909 keV) 
in both Chandra and XMM. Both the photon spectral in- 
dex and the normalization are fitted separately for each 
detector. Two absorption lines, one near z — and an- 
other near z = 0.03, were added following the above pre- 
scription. We required the parameters of the absorption 
lines to be the same in both detectors. Therefore, the 
model that we fitted to the data consists of a power-law 
continuum, two absorption lines, and (fixed) foreground 
Galactic absorption from cold gas. 

3.2. Continuum 

Upon fitting the composite model simultaneously to 
the XMM and Chandra spectra we obtain a spectral in- 
dex of 3.3 ± 1.7 and a 21.0-22.5 A flux of (7.76 ± 0.26) x 
10~ 13 erg cm -2 s _1 (unabsorbed) for the power-law con- 
tinuum of the XMM spectrum (90% errors quoted for 
both quantities). Because the Chandra LETG-HRC-S 
spectrum contains photons from all orders, not just the 
first order, some photons with wavelengths smaller than 
21.0 A in the higher-order spectrum will be mis-classified 
as having wavelengths between 21.0 — 22.5 A if the re- 
sponse matrix only includes information on the first- 
order diffraction spectrum. To account for the higher 
orders, the response matrix (<J2]) requires the spectral 
model to be defined down to wavelengths as small as 
- 5 A. 

We adopted the following procedure to define the con- 
tinuum for the Chandra data as a compromise between 
the desire to insure an accurate modeling of photons 
from higher orders and to measure the continuum us- 
ing only data very close to the energies of the absorp- 
tion lines. We initially fitted the power-law model mod- 
ified by foreground Galactic absorption from cold gas to 
the background-subtracted Chandra data over the broad 
wavelength range 5 — 30 A. This model provides a very 
good fit to the Chandra spectrum with a well-constrained 
power-law index, 1.883 ± 0.034. In all subsequent fits 
(though see to the narrow band (21.0-22.5 A) we 

restrict the the power-law index for the Chandra data 
to lie between these values established by the broad- 
band fit, thereby insuring an accurate modeling of the 
broad-band spectrum and thus the contribution from the 
higher orders. After fitting the composite model (i.e., 
continuum, absorption lines, foreground cold absorber) 
simultaneously to the XMM and Chandra spectra we 
obtain a 21.0-22.5 A flux of 5.52±° ;g x 10" 13 erg cm" 2 
s _1 (unabsorbed) for the Chandra power-law continuum; 
i.e, the continuum flux of the blazar measured by XMM 

10 http:/ /www. atm.ox.ac.uk/user/wells/voigt. html 



TABLE 1 
O VII Spectral Line Properties 







N 






Rcdshift 


(10 16 cm- 2 ) 


TO 


Milky Way 


-0.0004 - 0.0020 


(5.5,0.9) 


(20.3, 2.2) 


Sculptor 


0.0311 - 0.0327 


(4.7, 1.0) 


(18.7, 2.1) 



Note. — All parameters arc quoted at the 90% confidence level. The 
rcdshift ranges are determined by allowing the Doppler b-parameter to 
take any value between 1-300 km/s, a large range that should bracket 
physical WHIM values, and the redshifts themselves were allowed to 
take any values during the fits. The column densities and line-center 
optical depths arc lower limits evaluated assuming two values for the 
6-paramctcr (50, 100) km/s. 

is about 40% higher than observed by Chandra over the 
spectral region of study. 

3.3. Background 

Over 21.0 — 22.5 A the background comprises about 
8% of the observed XMM count rate and about 38% of 
the Chandra count rate. As noted above, we perform 
spectral fitting on the background-subtracted XMM and 
Chandra spectra of the blazar. However, below in ^3.41 
we also assess the statistical significance of the O VII ab- 
sorption lines using simulated spectra, for which we pre- 
fer to use model representations of both the source and 
background spectra for each satellite. We constructed 
models of the background spectra as follows. 

The Chandra LETG-HRC-S background spectrum 
over 5 — 30 A is fairly well described by a broken power- 
law folded through only the instrument RMF ("redis- 
tribution matrix file") and not the ARF ("auxiliary re- 
sponse file" ). This is achieved by defining a "BKNPOw/b" 
model in xspec vll. The model parameters are very 
tightly constrained. For reference, the best-fitting val- 
ues for the power-law indexes are, 1.552 and 1.888, and 
the best- fitting break energy is, 1.338 keV (9.27 A). 
Because the model is not a perfect fit to the data at 
all wavelengths, we adjusted the flux of the power-law 
to best match the spectral region of the oxygen lines. 
That is, adopting the values of the spectral indexes and 
break energy obtained from the broad-band fit we re- 
fitted the background model over the restricted range 
21.0 — 22.5 A of interest, which yields a normalization 
that is 4% ± 2% (la) less than obtained for the broad- 
band fit. 

For the XMM RGS1 the broad-band spectrum is not 
well described by simple models. Consequently, we 
adopted a single power-law model, again folded only 
through the RMF ("pow/b" model in xspec vll), but 
over a narrower wavelength region where the model is a 
good fit. We selected the wavelength range, 20 - 24 A, 
as the largest range producing a good fit while enclosing 
the critical 21.0 — 22.5 A interval. We obtain a power- 
law spectral index of 1.16 with a 90% confidence range 
(0.27-3.65). 

3.4. Absorption Lines 

In this section we discuss the significance and prop- 
erties of the absorption lines obtained from fitting the 
composite model to the data, while systematic errors in 
our analysis are examined in &]3.5l 

We list in Table [1] the redshifts of the two candidate 
O VII Ka absorption lines obtained from fitting the spec- 
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tral model simultaneously to the XMM and Chandra 
data. The line near 21.6 A is fully consistent with hot 
gas from the Milky Way (or Local Group WHIM), with 
the inferred error range on the redshift being consistent 
with zero. The redshift range of the line near 22.3 A 
matches very well the range of structure in the Sculptor 
Wall (z = 0.028-0.032) that intercepts the blazar's line- 
of-sight. Henceforth, we refer to these candidate lines as 
the Galactic and Sculptor WHIM lines. 

The change in the fit statistic, whether it be the C- 
statistic or x 2 , is very similar when adding the Galactic 
and Sculptor WHIM lines, indicating that the statisti- 
cal significance of the lines is also similar. For reference, 
the change in the x 2 statistic observed when adding the 
Sculptor WHIM line to the fit suggests the line is signif- 
icant at the rs 3cr level when interpreted in terms of the 
F-Test. However, the F -Test is not strictly ap plicable to 
spectral line detection (|Protassov et al.ll2002l ). 

To obtain a rigorous estimate of the statistical sig- 
nificance of the Sculptor WHIM line we employed a 
Monte Carlo procedure. Using as a reference the best- 
fitting spectral model excluding the Sculptor WHIM line, 
we generated with XSPEC synthetic XMM and Chandra 
spectra for the source and background ( £)3.3p with the 
same exposure times as the real data. The simulations 
were performed using the full energy resolution of the re- 
sponse matrices, after which the synthetic spectra were 
re-grouped as described in fJ2J In this way the simu- 
lations account for noise associated with the re-binning 
procedure. 

We fitted the synthetic background-subtracted spectra 
with the reference model and minimized the C-statistic to 
obtain a best fit. Then we added a second absorption line 
component with redshift range restricted to 0.028-0.032 
appropriate for the Sculptor Wall. After obtaining the 
new best fit, we recorded the reduction in the C-statistic. 
We performed this procedure for 10 simulations. The 
probability of a false detection is given by the number of 
simulations where the reduction in the C-statistic is at 
least as large as observed in the real fit relative to the 
total number of simulations. 

The fit of the reference model (i.e., without the Sculp- 
tor WHIM line) to the real data yields a best-fitting value 
of the C-statistic of 68.2 for a total of 88 spectral bins. 
Upon adding the Sculptor line and re-fitting, we obtain 
a C-statistic of 58.6 for a reduction of 9.6. For the Monte 
Carlo procedure described above, we find that in only 36 
out of the 10 4 simulations is the C-statistic reduced by 
at least 9.6. This translates to a false detection probabil- 
ity of 0.36%. Therefore, the statistical significance of the 
Sculptor WHIM line is 99.64%; i.e., about 3cr. Following 
the same procedure, but this time starting without the 
Milky Way line and with the Sculptor WHIM line al- 
ready in place, we determine the statistical significance 
of the Milky Way line to be 99.69%, also about 3o\ 

(We mention that the XMM data contribute more to 
the statistical significance than do the Chandra data. If 
the spectral fitting is performed separately for each data 
set, then following the above Monte Carlo procedure for 
the Sculptor WHIM line we obtain a statistical signifi- 
cance of 2.4ct for XMM and 1.7(7 for Chandra.) 

Apart from the redshift noted above, the line param- 
eters are not well constrained by the data. In Table Q] 
we list the 90% confidence lower limits obtained for the 



O VII Kq column density (N) and line-center optical 
depth (to = t{Eq,s)). Wc quote results for two values 
of the 6-parameter (50, 100) km s" 1 which span a large 
range of WHIM temperature and turbulent pressure. 

The lower limits of To ~ 1 for both lines demonstrate 
that the line centers are not optically thin. Therefore, 
one cannot infer the column density from the equivalent 
width independent of the fo-parameter; i.e., these lines do 
not lie on the linear portion of the curve of growth. For 
reference, the best-fitting equivalent widths are nearly 
30 mA for both the Galactic and Sculptor lines. 

The column density lower limits of w 10 16 cm -2 are 
consistent with expectations for both lines. Similar col- 
umn densities previously have been repor ted for the 
Mil ky Way (or Local Group WHIM) (e.g., iFang et all 
2003). As for the Sculptor line, we re-e xamined the 
cosmo logical hydrodynamical simulation of ICen fc Fand 
(2006) and identified five super-structures with similar 
(or higher) galaxy over-density relative to the Sculp- 
tor Wall. By shooting random sight-lines through these 
super-structures we obtain a mean O VII column density 
of » 3 x 10 16 cm~ 2 , consistent with the limits imposed 
by our spectral fits on the Sculptor Wall. 

3.5. Systematic Errors 

The large statistical errors obtained for the line col- 
umn densities and line-center optical depths dominate 
systematic errors associated with our data analysis and 
modeling choices, and therefore we do not provide a de- 
tailed systematic error budget for the line parameters. 
However, in this section we briefly discuss the possible 
impact of systematic errors on the estimation of the sig- 
nificance of the Sculptor Wall WHIM line, which is the 
principal result of our paper. We have in particular con- 
sidered the effects of errors in the background models 
and the wavelength calibration of the detectors, and we 
have also examined using \ 2 rather than the C-statistic. 

As we mentioned in ^3.31 the background level is only 
a small fraction of the total observed flux for the XMM 
RGS. Consequently, it is expected that small variations 
in the background level from that we adopted in our fits 
should have little effect on the significance estimate. In- 
deed, varying the background level by ±5%, which re- 
flects the statistical error on the power-law model ob- 
tained for the RGS background, did not have a noticeable 
effect. 

Because the Chandra LETG-HRC-S is unable to sep- 
arate spectra of different orders, we examined how much 
our results change if instead we used a response matrix 
containing only information on the first-order spectrum. 
Using only the first-order response leads to an underesti- 
mation of the actual background by almost 20% over the 
21.0 — 22.5 A range. Since the first-order response does 
not require a valid continuum model outside our range 
of interest, we allowed the slope of the continuum to be 
fitted without restriction over 21.0 - 22.5 A, which pro- 
vides a measure of the sensitivity of the results to the 
definition of the continuum. When performing simul- 
taneous RGS-LETG fits using the first-order response 
for the LETG we obtain a statistical significance of the 
Sculptor WHIM line of 99.53%, only slightly less than ob- 
tained when using the more accurate response containing 
orders 1 through 6. 

We also considered systematic uncertainties in the 
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wavelength calibration of the detectors. Current esti- 
mates indicate the calibrations are accurate to within 
2.4 mA in the RGS1 11 and to 10 mA in the LETG/HRC- 
S 12 . We find little difference in any of our results when 
allowing for shifts of these magnitudes between the de- 
tectors. We do mention that if the wavelength of the 
z = line is allowed to be fitted separately for the RGS1 
and the LETG then the best-fitting line centers are off- 
set by 48 mA. However, the shift is significant only at 
the ~ 90% confidence level, and it is substantially larger 
than the expected calibration uncertainty, indicating this 
marginal shift is probably noise. Therefore, we kept the 
line energies tied between the detectors. 

Finally, we also performed all of our fits by minimizing 
(data-weighted) x 2 an d obtained a consistent result for 
the significance of the Sculptor line as obtained using 
the C-statistic. In particular, the statistical significance 
of the Sculptor line obtained using \ 2 is, 99.81%, about 
3cr. 

4. CONCLUSIONS 

We report a detection at the 3<r level of an O VII 
Ka resonance absorption line located in the Sculptor 
Wall super-structure of galaxies. The column density (> 
10 16 cm" 2 ) of the line, though not precisely constrained, 
is consistent with that produced by similar structures 
generated in cosmological simulations. The presence of 
this absorption line is inferred from simultaneous analy- 
sis of XMM RGS1 and Chandra LETG/HRC-S spectra 
of the blazar H 2356-309 (z — 0.165). Since the line cen- 
ter is saturated (for both the Sculptor and Local lines), 
we interpreted the spectra with a physical model where 
the line is modified by the Voigt profile to account for 
both natural and Doppler broadening. 

The robustness of the detection benefits from three 
key aspects of our study. First, a critical factor in ob- 
taining this high detection significance level is that we 
know a priori the redshift of the intervening structure 
in the Sculptor Wall (z = 0.028 - 0.032). The impor- 

11 A. Pollock, XMM-Ne wton Users Group Meeting Pre- 
sentations, 6-7 May 2008, http://xmm.csac.csa.int/cxternal/ 
xmm_user .support / usersgroup / 200805 00/ 

12 M. Chung et al., http://cxc.harvard.edu/cal/Letg/Corrlam/ 



tance of this point was dramatized in the controversy 
over the claimed WHIM detection in the spectr um of 
Mrk421 (jNicastro et al.ll2005l ; iKaastra et al.ll2006l ). Sec- 
ond, the line is detected from a joint analysis of XMM 
and Chandra data. The consistent picture from the two 
satellites provides important reassurance given again the 
controversy for Mrk 421 where t he RGS did not confirm 
the claimed Chandra detection (Rasmussen "eTaII[2007h . 
Third, since the background source is a blazar, it is very 
unlikely that features intrinsic to the source contaminate 
our study of the foreground WHIM. 

This detection of WHIM in the Sculptor Wall rep- 
resents the most significant (non-local) WHIM detec- 
tion to date from X-ray absorption studies, providing 
vital evidence for the expected repository of the "miss- 
ing baryons" (e.g. , [Fukugit a et al.lll998t ICen &; Ostrikerl 
119991: iDave et alj |200lD complementary to that of X- 
ray emission studies o f WHIM in superclusters (e.g., 
IZappacosta et al. 2005) and the intersections of binary 
clusters (e.g.. IWerner et aLll2008[ L and O VI absorption 
line studies that probe WHIM at lower temperatures 
(e.g. JSavage et al.lll998l : iTripp et al.ll2006[) . 
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